On 
On 



Effective theory for the soft fluctuation modes 
in the spontaneously broken phase 
of the N-component scalar field theory 

A. Jakovac a 't A. Patkos 5 '} P. Petreczky 6 'Wd Zs. Szep 6 ' § 

(a) Department of Theoretical Physics, Technical University of Budapest, H-1521 Budapest, Hungary 
(b) Department of Atomic Physics, Eotvos University, H-1117 Budapest, Hungary 



O : February 1, 2008 

(N 
> . 

Q\ , Abstract 

m 

The effective dynamics of the low- frequency modes is derived for the O(N) symmetric scalar 
field theory in the broken symmetry phase. The effect of the high-frequency fluctuations is taken 
into account at one-loop level exactly. A new length scale is shown to govern the long-time 
asymptotics of the linear response function of the Goldstone modes. The large time asymptotic 
decay of an arbitrary fluctuation is determined in the linear regime. We propose a set of local 
Oh, equations for the numerical solution of the effective non-linear dynamics. The applicability of 

Q_i' the usual gradient expansion is carefully assessed. 

1 Introduction 

X 

In the study of dynamical phenomena with the participation of long-wavelength Higgs and 
Goldstone-modes one can account for the effect of the high frequency modes most conveniently 
by deriving a number of effective equations of motion. In a second step one might express the 
correction terms of the original equations with help of a set of auxiliary fields, allowing local 
representation of some nonlocal effects. A remarkable feature of this approach is that for pure 
gauge theories a simple kinetic interpretation can be given to the dynamics of the auxiliary 

fields |, !!!§. 

Recently we have shown that a similar interpretation is possible for fluctuations of the one- 
component self-interacting scalar field in the broken symmetry phase Q. The evolution of 
the modes with wave number |k| was considered under the assumption M < T; that is the 
mass scale set by the spontanous symmetry breakdown is smaller than the temperature. Since 
M ~ \/A( < ^), this condition relates the vacuum expectation value of the field to the temperature 
and is actually met in systems with second order phase transition not much below the transition 
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point. This condition actually was used when the comparison of the quantum results with the 
classical dynamics was made. 

The high frequency modes have a non-trivial impact on the dynamics of the low frequency 
modes if the scale characterizing the low-frequency fluctuations is much below the mass scale 
of the symmetry breaking. In this sense one has to go beyond the usual hard thermal loop 
approximation, since the effect of loops with momentum p ~ M < T should be taken into 
account. The classical statistical mechanical system proposed in || reproduces the linear source- 
amplitude response computed in quantum theory J?]]. 

In this paper we extend our discussion to the broken phase dynamics of the N-component 
scalar field theory with <I> 4 self-interaction. This model is relevant to the dynamics of the tt — a- 
system (N = 4) || and is actively investigated in connection with the phenomenon of the 
disoriented chiral condensate |], [R], [H], 12 1. The appearance, the evolution and the damping 



of the low frequency fluctuations and instabilities of the chiral order parameter are carefully 
studied in these papers. Less attention is paid to the damping of the Goldstone-modes. Though 
in the realistic case the explicit breaking of the 0(N) symmetry leads to massive "Goldstone" - 
bosons, it is of interest to see what is the intrinsic dynamics of these excitations in the "ideal" 
spontaneously broken case. It turns out that the damping of the Goldstone-fluctuations with 
frequency fco is essentially different in the respective domains 

M 2 

k « — « M < T (1) 

and 

M 2 

— « k « M < T. (2) 

Our conclusion is that both the on-shell dissipation time and the large time asymptotics of the 
Goldstone fields is qualitatively different in the first region in comparison to the fluctuations of 
the order parameter. 

The physical picture is very appealing. Radial fluctuations along the order parameter relax 
fairly quickly, what freezes fast the length of the vacuum expectation value of the <& field. The 
second stage of the relaxation consists of the slow rotation of the order parameter, which is 
described by the relaxation of collectively excited low momentum Goldstone fluctuation modes. 

The high-temperature one-loop quantum dynamics of scalar models has been investigated 
very actively recently partly as a kind of theoretical laboratory for developing powerful calcu- 
lational methods, partly with the aim to provide answers to important questions of inflatory 
cosmology and the physics of heavy ion collisions. Our work was substantially influenced by 
References ||, 0, |f, @, • 



In this paper we first derive in section 2 the effective equations for the low frequency modes. 
The coefficients of the terms correcting the classical equations are determined by various n-point 
functions of the fast modes. For the effective dynamics we work out the modification of the 
linear part of the equations (the self-energy operator), which is fully determined by the two- 
point function of the high-frequency fields. A detailed study of this quantity appears in section 
3. Stated in more technical terms, we compute one-loop contributions to the two-point functions 
with full T = propagators and restrict the temperature dependent contribution to the po > A 
modes. This approach follows the finite temperature renormalization group transformation 
scheme of D'Attanasio and Pietroni |l8|]. In section 4 we shall discuss, in particular, whether the 
linear response of the high frequency modes to long wavelength fluctuations allows a classical 
kinetic theory interpretation. Using the results obtained for the two-point function of the 
theory, in section 5 we present the explicit effective field equations, with help of auxiliary fields 
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introduced to handle the nonlocal nature of the effective dynamics. Our results are summarized 
in section 6. In Appendix A some results of the main text are rederived with help of the 
conventional perturbation theory in a form explicitly showing its equivalence to the Dyson- 
Schwinger treatment. Appendix B presents the iterative solution of the dynamics of the classical 
0(A)-model. Some relevant integrals are explicitly evaluated in Appendix C. 



2 The effective equations of motion for the slow modes 



The Lagrangian of the system is given by 



L 



1 



1 



24 



((<Pc 



\2\2 



(3) 



The aim of this investigation is to integrate out the effect of po > A high frequency fluctuations. 
From the point of view of the final result it turns out to be important whether A < M or A > M, 
where M is the mass scale spontaneously generated in the broken symmetry phase. In the spirit 
of the renormalization group the value of A will be lowered gradually and the importance of 
passing the scale M will become evident in this process. Our actual interest will concentrate 
on the temporal variation of the lowest frequency fluctuations ko « A. 

We separate in the starting Lagrangian the high-frequency modes (4> a (x) with frequency 
<jj > A) and the low-frequency modes (<& a (x) with frequency u < A) 

Mx) -> $a(x) +M X )- ( 4 ) 

Averaging over the high frequency fluctuations gives (<j) a (x)) = which yields also & a (x) = 
((p a (x)). This separation leads to 



L = \(d*a) 2 + \{d<t) a ?- r ^[(S>a) 2 + 



24 



{$>a?{$>b? + 



24 



2 



(5) 



4($ O O )($ 6 6 ) + 4(^ a <P a )(<P b y + A(^ a 4> a )(1> b y + 2($ a )^ fe ) 



From Eq.(||) the equations for the slow modes can be derived: 
(8 2 +m 2 )4> a + A [8M^b) + 4$ a ($ b ) 2 + 40 a (^) 2 +4(/)a(^) 2 + 8$ a (^) + 4$ a ( 



0. (6) 



The effect of the high-frequency fluctuations on the slow ones is obtained by averaging the 
equations with respect to their statistics. At one-loop level accuracy (0 a 060c) = and one 
arrives at 

52 , , ^SL r^_\2 , I ± i \ , It ±_ \2\ 



(d l + m^ a + -$ a ($ 6 r + ^b(Mb) + ^a((<f>b) 
DO D 

We introduce at this point the conventional notation: 

A a bO,y) = {<p(x)<p(y)). 



0. 



(7) 



(8) 



Below no summation will be understood when repeated indices appear without the explicit sum- 
mation symbol. We shall also use specific pieces extracted from the above two-point functions 
defined as 



A i?( x >y) = A ab(x,y)\^ =0 , 



(9) 
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In the broken symmetry phase one has to separate the nonzero average value from the slowly 
varying field, and write the equation only for the fluctuating part: 



*a-»**al + *a, (10) 

where we have chosen the direction of the average to point along the a = 1 direction. We shall 
analyze the resulting equations for the a = 1 and the remaining a = i ^ 1 components separately. 
Since the main effect of the high frequency modes we are interested in is the modification of 
the mass term (self-energy contribution), therefore we shall restrict our study to the linearized 
equation of Q a (x). One has to take into account that the two-point functions {(j) a <t>b) depend on 
the background <& a {x)- For the linearized equations it is sufficient to compute them only up to 
terms linear in the background. 
The two equations are 

(a 2 + \^)^{x) + ^\{x) + *< (*) + E + J i (*) = o, 

dH i (x) + ^ i {x)^ l (x) + ^ i {x)Y j ^l{x) + Ji( X ) = 0, (i^l) (11) 

with the induced currents 

Ji(x) = ^$A$(x,x) + ^(N - l)*A$(x,x), 

J 4 (x) = ^(Ai 0) (x,x) - Ag>(*,*))*, + ±*A$(z,x). (12) 
We have simplified the equations using the fact Aj^j has no ^-independent piece, since with 



zero background the correlators are diagonal. The equations (12) express the linear response of 
the fast modes to the presence of a low frequency background producing an effective source term 
to their classical dynamical equations. Also one has not to forget that $ is now the solution of 
the constant part of the equations: 



m 2 $ + ^(3Aff(x,x) + (N — l)A$\x,x))Q + ^<S 3 = 0. (13) 

3 The two-point function of the fast modes 

Varying Eq.@ with respect to <f)(x) one arrives at the equations of motion of the fast modes in 
the background of <&(x): 

(d 2 +m 2 )<£ a +A(8$ a ($ 6( £ 6 ) + 4^^ = 0. (14) 

The equations of motion linearized in the high frequency fields can be written in the form 

j(<9 2 + m 2 )5 ab + ^[2$ a $ b + ($ c ) 2 6 ab }\ct> b = -^ a ($ b ) 2 . (15) 

We introduce the notation 

m 2 ab (x) = m 2 5 ab + ±[2* a (x)* b {x) + {$ c (x)) 2 5 ab ] (16) 
6 
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and apply Eqs.(15) and ( |l6| ) to the fields appearing in the definition of the two-point function 

A ac (x,y): 

dlA ac (x,y) = -m 2 ab (x)A bc (x,y), 

d 2 A ac (x, y) = -m 2 cb (y)A ab (x, y). (17) 

In the derivation of these homogeneous equations we have exploited that (<b a (x)(<b c (x)) 2 (pb(y)) 
= 0. After the Wigner-transformation 

A(X,p) = J d\e ipu A{x,y), u = x-y (18) 

one arrives at the exact linear equations for the Wigner transforms. (The use of the Latin 
letters k,p,q is reserved for the 4-momenta and pk,kq, etc. denotes their Minkovskian scalar 
products.) For the two-point functions diagonal in 0(N) indices one finds in the broken phase 
to linear order in the background (after using Eq.([lO|)) 

[\d 2 x -ip-d x -p 2 + M 2 a ]A aa (X,p) + \ a § J A*Q 1 ( q )A aa (X,p-l)e- i * x = 0, 
l\d 2 x + i p-d x -p 2 + M 2 ]A aa (X,p) + \ a $ J ^ i $ 1 ( 9 )A oa (X,p+|)e-^ x =0, (19) 
where for the Goldstone bosons on has 

M 2 = 0, Xi = ^, (20) 

while for the heavy ("Higgs") mode 

M 2 = h 2 , Ax = A. (21) 

Here the tree level mass relations were used, as they correspond to the actual order of the 
perturbation theory. This procedure can be improved using Eq.(13) for what is equivalent 
to the resummation of the perturbative series. 

After performing the Fourier-transformation also with respect to the center-of-mass coordi- 
nate and subtracting the resulting two equations one arrives at 

2pkA 11 (k,p) = -A# J -|$i( g )[A u (A - q,p + |) - Au(fc -q,p- |)], 

2pkA u (k,p) = ~^J ^^i(q)[Aii(k - q,p + |) - A u (k - q,p - |)]. (22) 

For the mixed two-point function An in the same steps a similar equation is derived: 

(2pk + ~$ 2 )Aa(M =-~*J 704 - q,P + |) - An (A; - g,p - |)]^(fe). (23) 

In the kinematical region, where the center-of-mass variation is very slow one can assume 
q~k. (If there would be no X-dependence one would find A{k — q, P) ~ 5{k — q)). If in addition 
one restricts the variation of p in the high frequency fields by the relation | fc„ | < < A ~ |p„ | , 
one can approximate the integrals in these equations by a factorized form. If the functions 
A aa {k,p — |) are replaced by the first two terms of their power series with respect to k/2, one 
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recovers the drift equations proposed originally by Mrowczynski and Danielewicz |]13| . In case 
of the diagonal two-point functions these equations are identical to the collisionless Boltzmann 
equation for a gas of scalar particles, whose masses (M 2 + \<&&i(x) and Mf + A$<I>i(x)/3) are 
determined by the <3?i field: 



P ■ d X A aa (X,p) + ^Ox^dpAaaiX^) = 0. 



(24) 



Below we shall discuss an alternative to this expansion, which does not exploit the assumption 
of slow X-variation. At the end we will be in position to assess the validity of the assumption 
which led to the Boltzmannian kinetic equations (|24|), 



At weak coupling one can attempt the recursive solution of Eqs. (22) and (p23[). The starting 
point of the iteration is the assumption that the high frequency modes are close to thermal 
equilibrium, therefore one has for the starting 2-point functions 

A^(p) = 2n5(p 2 )(e(p )+h(\p \)), 

Ag>(p) = 2nS(p 2 - M 2 )(0(p o ) + n(bo|)), 

1 



n(x) = n(x)Q(x — A) 



3 /3x _ i 



<d(x - A). 



(25) 



Since the starting distributions are independent of X, the integrals in Eq. (19) factorize exactly 
and one has for the first corrections the following explicit expressions: 



= -^A.**i(*)Kg?(p + ±) - AS(p- -)], 



AS } (fc,p) 



a|?(p-5)]. 



(26) 



The retarded nature (ie. forward time evolution) will be taken into account with the Landau 
prescription ko — ► ko + ie (see Appendix [A| for the equivalent formulas in the conventional 
perturbation theory). 



4 Non-equilibrium linear dynamics of the slow modes 



The quantum improved, linearized equations of motion are (see eq. (|TT 

(8 2 + M 2 )<Z> a (x) + J a (x) = 0. 
The (retarded) self-energy function is introduced through the relation 

J a (x) = Jd 4 yU a (x-y) $ a (y), 

where, after Fourier-transformation one obtains, (see eq. (|2 

d 4 p 

pk 



ni(fc) 



Ui(k) 



d A p 

WT 4 



+ 



A$ 



(jv - 1)^ (aJ?) (p _ |) _ a2> (p + 



A$ 



1 



A (0 V 

3 / 2^ + M 2 ^ n[P 



4(Ag)(p)-Aif)(p) 



2 ; ii ^ T 2 



(27) 
(28) 



(29) 
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These integrals can be calculated (or at least reduced to ID integrals). Some details of the 
calculations can be found in Appendix [C], here we just state the results: 



III (A;) 



^f) Rx^MJ + t^-) (N- 1)^,0) 



Mk) = ( —3— 1 Ri(k,Mt] 



(30) 



The explicit form of their real parts is the following: 

00 

ReR 1 (k,M) = Jdu) (1 + 2h(u)) 

(l + n{u))A( 
(l + n(|p|)M( 



'a-^ tj ) + a( 2|p||k 



CO 



ReRi{k,M) = J duj 

M 



2cok + k 2 ■ 
2|p||k| 



+ /d|p| 



2ujk + k 2 + M 2 
2|p||k| 



-2u;k + k 2 
+ n(u))A{ 



2|p||k| 



-2|p|fc + k 2 - M 2 



)+n(\v\)A( 



-2ujk + k 2 + M 2 

gjpjjkj 
2|p|fc + k 2 - M 2 ' 



(31) 



where |p| 2 = u> 2 — M 2 and 



A(x) 



1 



-arth(x) 



1 



4vr 2 |k| v ' 8vr 2 |k 
The expressions of the imaginary parts look as follows: 

-1 



In 



1 + x 



1 — X 



(32) 



InxRi(A;,M) 
ImRi(k,M) 



Pc r. 



47r|k| 

-1 
16vr|k| 



Q(-k 2 )J dPh(P) + ®{k 2 -AM 2 ) J dP(l + h(ko-P)+n( p )) 

Pc-ko k /2 

Q++fco |<3- 1 

0(M 2 -k 2 ) J dPh(P) + J dPh(P) 

Q+ \Q--ko\ 



-@(k 2 - M 2 ) J dP (1 + h{k ~P)+ n{P)) 



where 



|k| 
~2~ 



1 



4M 2 & 
1 — - 

fc 2 2 ' 



fc 2 -M 2 



2(A: + |k| 



Q- 



k 2 -M 2 
2(k - |k 



(33) 



(34) 



The real part is divergent at zero temperature which cancels against the coupling constant 
counterterm contribution. Care has to be taken, however, when implementing the regularization, 
to maintain the Lorentz invariance for the pieces containing T = parts of the propagators A^°\ 
which is manifest in the original form in Eq. (pS|). 

It is remarkable that the usual domain of Landau damping (|k| 2 > \ko\ 2 ) is apparently 
extended up to M 2 + |k| 2 > |/cq| 2 - An analogous situation has been noted and interpreted 



recently for the propagation of a light fermion in heavy scalar plasma [20]. Below we find for 
the damping of soft on-shell Goldstone-modes the same interpretation. 
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The linearized equations of motion can be analyzed from several points of view. One is 
the determination of the dispersion relations. This describes the quasiparticles, and physically 
corresponds to the "dressing" of an (external) particle passing through the thermal medium. 
The other point of view is the field evolution when one follows the solution of the field 
equations developing from given initial conditions (history). 



4.1 Dispersion relations for on-shell waves 

The position of the poles is determined by the equation 



k 2 -M 2 a -Yl a {k) = 0. 



(35) 



We split ko into real and imaginary parts: k$ 
perturbative solution can be written as 



uJ 



ujq + ReII(u;o, k) 



iT, and assume T -C uj. Then the 
Indigo, k) 



2lj 



(36) 



where u>q = M 2 + k 2 . The corresponding time dependence for fixed wave vector k and given 
initial amplitudes (d t $(t = 0, k) = P(k), $(f = 0, k) = F(k)) is found 



*(t,k) 



, , smut , . 
P(k) + F(k) cos cot 



-rt 



(37) 



Goldstone modes. The tree level dispersion relation has no mass gap. The second equation 
of Eq. ( p9[ ) shows that IT (A; = 0) = 0, that is no mass gap is created radiatively neither at 
finite temperature; this is the manifestation of the Goldstone theorem. The on-shell imaginary 
part, as Eq. ( |30[ ) shows, comes from the continuation of the Landau-damping extending up to 
k 2 < M 2 . Going back to Appendix C (Eq.( |lli| )), one finds that the imaginary part receives 
contribution from the collision of a hard thermal Goldstone particle with distribution n(po) 
with the soft external Goldstone wave of momentum k producing a hard Higgs-particle minus 
the inverse reaction, when a hard thermal Higgs with momentum distribution n(po + fco) decays 
into a soft and a hard Goldstone. The two contributions can be combined into a single integral 
leading to the following expression for the damping rate: 



r<(k) = — 



i 



3~ ) 32vrk 2 



M?/4|k|+|k| 



dph{p). 



(38) 



Aff/4|k| 



If |k| <C M\ we can write with a good approximation 




32vr|k| l 4|k| ; ' 

This contribution survives the IR cut, if A < M 2 /4|k|. 

The result (^) is interesting from several points of view, 
where we neglect all the masses, we would found Tj ~ ^(0) 
classical approximation corresponds to the substitution n(x) —■ 
(see Appendix |B|) in 



r 



class 



AT 
247' 



(39) 



First, in HTL approximation, 
= 0. On the other hand the 
T/x in Eq . (j39|) , which results 



(40) 
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However, for very small momenta (|k| << M 2 /4T) the correct result is exponentially small, de- 
viating considerably from the classical result. The purely classical simulations therefore cannot 
reproduce the correct Goldstone dynamics in the long wavelength region, they would overesti- 
mate the damping rate. 

The most important consequence of the form of the Goldstone damping is the generation 
of a new dynamical length scale M^ 1 = 4T/Mf . The components of the Goldstone condensate 
with longer wave-length can not (exponentially) decay, they survive for a longer time. 

Higgs modes. From Eq. (|29|) one can easily read off the on-shell (k^ = k 2 + M\ ) value of 
the imaginary part of the self-energy. This is entirely due to Higgs scattering into a soft + hard 
Goldstone-pair, and leads to 

|k|/2 

-Imni AM, 2 ,„ _ „,k n _ ,kn 



2kQ 967r/co|k 
In the k <C M\ limit it simplifies to 



±- (N - 1) J dP (1 + n(| - P) + n(y + P)). (41) 



r 



(N-l)(l+hA). (42) 



AMi 

96vr v " ' V2 2 



This survives the cut if 2A < Mi . 

Since M\ < T we can perform a high temperature expansion. This means at the same time 
that the classical approximation is applicable. We find 



AT 

r^-fA-D, («) 



which is (N- 1)/2T; 



class 



4.2 Solution of the initial value problem of the slow fields 

Equation(^Tj) describes an integro-differential equation. For its solution we have to know the 
complete past history of the fields. The previous subsection provides a specific way to look at 
the problem. In general we can introduce z a {x) = J yo<0 d 4 yH a (x — y) & a (y), then for xq > 

Since all singularities are in the lower half plane, & a (%) defined by this relation vanishes for 
xq = t < 0. z{x) can be used for setting the initial conditions |Q, as it can describe jumps in 
the field as well as in its time derivative. For example & a (x) = for x$ < and $ a (t = 0) = F a , 
dt$ a (t = 0) = P a corresponds to 

z a {x) = -P a (x)6(x ) - F a (x)5'(x ), z a (k) = ik F a (k) - P a (k). (45) 



Direct substitution of Eq. ( f44| ) into Eq. ( p7| ) shows that it satisfies the effective homogenous 
wave equation. By evaluating the ko integral for t = (for details, see below) one also can 
demonstrate that it fulfills the initial conditions set above. 

In general, we are faced with the computation of Fourier transforms of functions analytic 
on the upper complex plane. The procedure of extracting their large time asymptotics was 
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investigated already in Refs. [Q, For t > the k$ integration contour can be closed with an 
infinite semi-circle in the lower half-plane and we pick up the contribution of the cuts and poles 
inside the closed integration path: 



fit) = j ^ f(ko) = £ + E /|| P(ko) e~^, (46) 

— OO U) epoles /ecuts j 

where Z(u) is the residuum of the physical poles of the /-function [], p = z'Disc/ is the discon- 
tinuity along the cut, / = [wi,^] is the support of the cut. The second term can be computed 
again by completing the original integration interval to a closed contour. The discontinuity 
itself may have poles (but no cuts!) which contribute in the same way as the "normal" poles. 
After these poles are "encircled", there are two straight contours - parallel to the imaginary 
axis - left, starting at the two ends of the cut and running in the interval [u,uj — zoo], where u 
is the end (or starting) point (threshold). After this analysis one arrives at the generic form 

oo 

f(t) = £ (-iZ(u,))e+* + £ {T)e^ j^p{u-iy)e-y\ (47) 
^epoles ^ethresh. o 

where the — sign is to be applied for the starting, the + sign for the end point of the cut, and we 
have to take into account all poles, also the ones on the unphysical Riemann sheet. Expanding 
p around u into power series of y, the y integration can be performed. The term ~ y n of the 
expansion contributes to the time dependence t~ n ~ 1 . The large time behavior is dominated by 
the lowest power term of the expansion. If p cannot be power expanded then the damping for 
large times is faster than a power-law. 

In our case the position of the poles is determined by the dispersion relations (see previous 
section), and for given initial values we find 



%ole{t,k) =Z(k) 



„ „ . sinc^i . 
P(k) +F(k) coswt 



e- Lt . (48) 



This expression coincides with eq. ( |37| ) apart from the wave function renormalization Z(k). The 
time dependence of the quasi-particle pole contribution was analyzed before. 

For the cuts the factorized form can be used for the spectral function: p a = z a Pa, where p„ 
is the spectral function of the propagator: 

G _ -2Imn a 
Pa ~ (k 2 — M 2 — Ren a ) 2 + (Imn a ) 2 ' ( 9) 



The thresholds and the leading large-time behavior are determined by Imll, Eq. (p3[). The 
imaginary part of the Goldstone self energy is non-analytic only at k 2 = 0, but also there the 
non-analytic piece vanishes as ~ exp(— M 2 /(2T\ k$ — |k||)). This finally leads to a damping 
faster than any power which can be seen after a saddle point analysis 



1/4 



cut A 



XM 2 \ t 3 



, sin I kit _ ., . 
P(k) .J , 1 + F(k) cos k t 



(50) 



^^The poles below a cut do not contribute to this formula, they lie on the unphysical Riemann sheet. 
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The cut contribution for the Higgs damping, on the other hand, is similar to the zero temperature 
case 



*«rt,l(*» k ) = - 



AT 1 
24Mi (ttM^) 3 / 2 
A|k| 1 



P(k) 



sin(2Mit - it/4) 



96vrAfi TiMit 



P(k) 



2Mi 
sin(|k|t + tt/2) 



+ F(k) cos(2Mit - vr/4) 



+ F(k) cos(|k|i + vr/2) 



(51) 



The first term comes from the threshold of the Higgs pair production, the second from the 
threshold of the Goldstone pair production or their Landau damping. The Landau damping of 
the Higgs particles does not contribute to the power law decay. 

The terms decaying as some powers of time will dominate the time evolution after the period 
of the exponential decay for the Higgs bosons (|k| _1 <C M^ 1 ). Similar is the case for small 
Goldstone domains (|k| _1 <C M^ 1 ), there the exponential decay is followed by a ~ exp(— ^/t) 
time evolution. In case of large Goldstone domains (|k|~ x 3> M^ 1 ), however, because of the 
exponentially small damping rate, the situation is reversed: the ~ exp(— y/t) behavior will be 
dominant for intermediate times, while the amplitude is reduced by a factor of Z . Only for 
very long times will become the exponential damping term, arising from the Goldstone-pole, 
relevant. Its action will erase completely the large size domains. 



5 Nonlinear dynamics 

The calculation of the effective non-linear evolution of the low-frequency modes can be based on 
Eq. ( |Tl| ) using the explicit expressions for the Fourier transforms of the induced currents defined 
in Eq.(|l|): 

+ (Ar_1) ^? / ^^iW[4 0) (p-fc/2)-Al 0) (p + fc/2)], 

J + k) = ~-3 I ^^W^^ ~ k,2) ~ A ^ + (52) 

Implicitly in all Wigner-transforms A^ ) the Bose-Einstein factor is understood with a low 
frequency cutoff, well separating the modes treated classically from the almost thermalised high 
frequency part of the fluctuation spectra. 

For the purpose of numerical investigations the above non-local form of the induced currents 
is difficult to use. In this section we discuss the introduction of auxiliary fields making the 
numerical solution of the nonlinear dynamics easier to implement. We shall not take into 
account the nonlinear piece of the source induced at one-loop level, since in weak coupling the 
leading nonlinear effect comes from the tree-level cubic term. The consistent inclusion of the 
higher power induced sources will not lead to qualitatively new, leading effects unlike the linear 
source, which is responsible for damping. We leave this extension for future investigations. 

We concentrate first on the induced current J\. By appropriate shifts of the integration 
variable p it can be rewritten as 



11 



After explicitly performing the po integration one arrives at the following expression which has 
a clear interpretation specific statistical average: 

-A 2< l> 2 f d 3 p 



Ji(fc) = $i(k) 



1 



+(N-1) 



A 2 <P 



(2tt)3 2p 
d 3 p 1 



l + 2n(po))( 
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■(l + 2n(po))( 



1 

2p£; + k 2 
1 



1 



2^fc2 Jl w=(P 2 + M i 2 ) 1/2 



1 



(54) 



(27r) 3 2po Vi ' "'^""^pk + k 2 2pk-k 2,lpo=M y 

The temperature independent piece corresponds to the T = renormalization of A, fixed at the 
scale k 2 . It is absorbed into the mass term of <&i, therefore we retain in the integrand only the 
terms proportional to the cutoff Bose-Einstein factors. 

We introduce two complex auxiliary "on-shell" fields W a (x, p), a = They correspond to 
the two different mass-shell conditions appearing in the above equation, and fulfill the equations: 



(2pk- k 2 )W a (k,p) = $i(fc). 



(55) 



Clearly, J\ (x) can be expressed as a well defined combination of the thermal averages of these 
fields: 

N 



Jix) 



A 2< i> 2 



((W l (x,p)) + (W\x,p) 



1 



9 



((W\x,p)) + (W\x, P y))). (56) 



Thermal averages are defined by the usual formula 

d 3 p 1 



(W a (x,p)} 



(2vr) 3 2p 



n(Po)W a (x,p)\ po={p2+M 2 ) i/2, 



hi. 



(57) 



For the relevant combination one can use in place of Eq. (55) the equations arising after the 
combination W a (x,p) — W a (x,p)* is eliminated: 



[(2 P ky - (k 2 y](w a (k,p) + w a (-k, P )*) = 2k l $ 1 (h). 



(58) 



Equations (11), ( pq ) and ( |58| ) represent that form of the nonlinear dynamics which is best 
adapted for numerical solution. 

For very small values of k$ a simplified form can be derived which coincides with the result 
of the conventional kinetic treatment of the Higgs-modes ||. One arrives at this approximate 
expression of J\{k) if one performs an appropriate three-momentum shift p — > p =F k/2 on the 
p variable in Eq. (|54|). One finds the following expressions for the denominators, accurate to 
linear order in ko: 



2pk ± k 2 -> 2A:ov / (P : Fk/2) 2 + M 2 - 2pk ± k^ » 2pfc(l ± k /2p ) + 0{ 



"0 ' K / > 



Po Po =F 



pk 

2po' 



(59) 



Performing the expansion of the denominators and of the arguments of n(po)/po to linear order 
in k, &o ° ne finds the following approximate expression for J\: 



Ji(k) 



I 



d 3 p 



dpo pk 

+ if I 

pk ' po V 2 



+ ^(l+"(P0 



■^0 dp • 



jV-1 i / dn(pp) kopp , J_/l , ^f., \ 



dn(po) 



lpo=|p| 



po=(p 2 +M 1 2 ) 1 /2 



(60) 



The IR-cutoff A >> |k| is important to prevent the singularity of the contribution from the 
Goldstone- modes. The contributions from the second terms in each line are independent of 



2 The proper solution of the effective dynamics of the heavy field surrounded by massless Goldstone-quanta will 
require the application of methods analogous to the Bloch-Nordsieck resummation at T = We thank M. 

Simionato for an enlightening discussion on this point. 
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k, therefore they are simply understood as a finite shift in the squared mass of <3?i. With these 
steps one arrives at the final expression for the nonlocal part of the induced "Higgs" current: 



Ji(k) 



d?p r 1 dn(po) 

W?^pI~^po~ pq ~- 



N — 1 1 dn( Po ) 



1 

vk 



k ^(k). (61) 



Here = (l,p/po)- In this case it is sufficient to introduce only two real auxiliary fields 
W 1 (x,v) and W l (x,v) in order to make the dynamical equations local. They are defined 
through the equations: 

vkW a (k, v) = k $i(k), a = l,i. (62) 

The weighting factors in the integrals over the auxiliary field can be interpreted as the deviations 
from equilibrium of the Higgs and Goldstone particle distributions, <5ni(p) and <5nj(p), due the 
scattering of the high frequency particles off the $i condensate. Then one writes 



J l (x) = - — (5(W 1 (x,v)) 



+ ^_l S (W i (x,v)) 



(63) 



where §{■■■} means a phase space integral weighted with the nonequilibrium part of the relevant 
distributions. 

For this case it is easy to construct the energy-momentum vector corresponding to the 
nonlocal piece of the induced current. One can follow the procedure proposed by Blaizot and 
Iancu [22] and investigate the divergence of the (/x, 0) component of the energy-momentum 
tensor: 

dA ced = -Jld^i, (64) 

The task is to transform the right hand side into the form of a divergence. Using Eqs.(|6l|) and 
j2D one finds the following expression: 



-Ji9°$i 



+ 



J ^(w\ X ,v)j:^{vd x )W\ X ,v)\ 
^W\x,v)y^(vd x )W\x,v)\ po= ^ 



po = (p2 + Af 1 2 )l/2 



(65) 



From here one can read off the induced energy-momentum function of the approximate Higgs 
dynamics: 



rp^.0 

induced 



A 2 <£ 2 f d 3 p r 1 dn(p ) 



N 



2vr) 3 
1 1 dh(pQ 



Po d Po 



p\ dPo 
-VH*,v) 2 | p(Hp 



Wl (^ V ) 2 U = (p 2 +Af2)l/2 



(66) 



For the exact (essentially non-local) dynamics we did not attempt the construction of an 
energy functional, but its existence for the above limiting case hints for the Hamiltonian nature 
also of the full one-loop dynamics as expressed in terms of the auxiliary variables. 

The analysis of the induced Goldstone source, appearing in Eq.(|52|) goes in fully analogous 
steps. The temperature dependent part to which a statistical interpretation can be linked is 
rewritten with help of two complex auxiliary fields V a (x, p), a = l,i as 



Ji(x) 



d z p r 1 



(2vr) J 



-n(po = (p 2 + M 2 ) 1 / 2 )R e y 1 (x,p) - — n(po = |p|)Re V\x, p) 
Po Po 



(67) 
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where the auxiliary fields fulfill the equations 

(2pk + k 2 + Ml)V l {k,v) = (2pk + k 2 )^ l (k), 

(2pk — k 2 + M 2 )V i (k, p) = (2pk - k 2 )<5>i(k). (68) 

In this case even for very small values of k$ we were not able to derive any limiting case in 
which one could treat the time evolution of the low frequency Goldstone modes as a truly 
Boltzmannian kinetic evolution. 

The solution of the system of equations ([ll]) , (pq) , (|58|) , (|67|) , (|68|) requires the specification 
of initial conditions also for the auxiliary fields W a + W a * , V a . If one uses the "past history" 
condition for the physical fields $>i(x) = const, <3?i(x) = const for t < 0, then the linear integral 
equation form of the equations for the auxiliary fields and their retarded nature imply vanishing 
W a , V a for t = 0. 



6 Conclusions 

In this paper we have performed a complete 1-loop analysis of the time evolution of low frequency 
(fco << Mi) field configurations in the broken phase of the O(N) symmetric scalar field theory. 
We have shown explicitly the full formal equivalence of the leading order iterative solution of the 
Dyson- Schwinger equations and the perturbation theory computation of the two-point function. 

We have analyzed explicitly the case when T > Mi, which occurs, for instance, in the 
vicinity of the second order phase transition of the model. Here we could compare our results 
with the results arising from the dynamical equations of the classical O(N) field theory with 
appropriately chosen parameters. For the "Higgs" particle we have found complete agreement 
for the |k| << M\ < T modes. Related to this is the fact that we were able to show that 
the exact one-loop equations derived for the auxiliary fields W a have a simple Boltzmannian 
collisionless kinetic form for small |k|. 

However, the analysis gave different conclusions for the linear response of the Goldstone- 
modes. The agreement with the classical theory is restricted to the interval M 2 /4T « |k| << 
M±. Below the new characteristic scale 

M 2 j± (69) 

the damping of the Goldstone modes becomes exponentially small and it vanishes non-analytically 
for |k| — > 0. This is a very natural manifestation of the Goldstone-theorem in a dynamical sit- 
uation: no homogeneous ground state will relax to any "rotated" nearby configuration. In the 
light of this suggestive picture it is not surprising that we could not find a classical kinetic in- 
terpretation of the exact one-loop dynamics of the Goldstone modes even for very small values 
of |k|. 

Besides the exponential damping analyzed above, (1 — Z) fraction of the initial configuration 
follows a different time evolution determined by the particle production and Landau-damping 
cuts. In case of the lowest wave number Higgs fluctuations this leads to the result that the 
exponential regime will be followed by a power-decay for times t > M-j -1 . In case of the 
Goldstone modes the cut contribution is ~ exp (—Ay/M^i), and for modes above the scale M 2 
this will dominate for large times. Below the scale M 2 the cut contribution will be observable 
for intermediate times, while the pole dominated damping, because of the exponentially small 
damping rate, becomes relevant only for very long times. 
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It will be interesting to see through numerical investigations, if the onset of the nonlinear 
regime will influence the damping scenario of the Goldstone fluctuations. The effect of interac- 
tion among the high-A; modes (in addition to their scattering off the low-k background) on the 
effective theory merits also further study. Finally, we work on the extension of our analysis to 
the Gauge+Higgs models in the broken phase, relevant to the physics of the standard model 
below the electroweak phase transition. 



A Appendix 



Instead of the generalized Boltzmann-equations described in Section |3| we can use also pertur- 
bation theory to evaluate the two-point correlation functions. Here we will demonstrate the 
equivalence of the one-loop perturbation theory and the iterative solution of the generalized 
Boltzmann-equations in the case of A^ 1 ). 
In the perturbation theory we write 

<<^)<M*)) = i (T^^tf^ 



-iSi 



where </> ( - ' are the free fields, Sj is the part of the action which decribes the interaction between 
the different field components in the presence of the background. T c stands for the time ordering 



along a complex time path c specified in [19]. At one loop, to linear order in <3? we need only 



Si 



IT 



Oc 



d 6 x 



$1 {4> b f + 2($ 6 6 )0i 



where the integration variable xq c represents the points on the complex integration contour in 
the t plane. The field operator contractions are performed with help of the matrix propagators 



iG> b (x) 



iG< b (x) 



<T0 o (x)0 6 (O)> (M0)Mx)) 



(70) 



where T* denotes anti time ordering. The tree level propagators are diagonal G ab (x) = 5 ab G a (x). 
Since the background depends on the real (not the contour) time we find 



(Mx)Mx)) {l) 

where 



x 



do) 



(x)Sj 



d 4 y [$i(y)6 ab + <f> a (y)5 b i + <f> b (y)5 al ]S ab (x-y), 

(71) 

(72) 



iS ab (z) = G%(z)GC(z) - G<(z)G<(z). 
The propagators can be expressed with help of the spectral function p as 

G<(p) = n(po)p(p), G c (t,p) = e(t)p(t,p) + G<(t,p). (73) 
Using the (t, p) representation and finally performing time Fourier transfromation leads to 

d 3 p dp dp' p a (p , p)p b (p' , p + k 



Sab(k) 



(27I") 3 2-7T 2"7T ko — P0 — Po + ^ 



(l+n(po) +n(Po))- 



(74) 



Using free spectral functions p a (p) = (2Tr)e(po)5(p 2 
(see for example ]TiJ|) 

d 3 p 



we arrive at the known expression 



Sab{k) 



1 

(2tt) 3 4u; a u b 
n a - n b 



l + n a + n b 



k 
+ 



1 + n a + n b 

uj a - uj b + ie k + uj a + u) b + ie 
n b - n a 



+ 



k + io a — io b + ie k — co a + u> b + ie 



(75) 
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On the other hand introducing in (|74| ) A^(p) = G > (p) = (1 + n(p ))p(p) = G < {—p) 1 next 
performing one of the po integrals, and shifting properly the p integral we obtain 



Sab(k) 



d 4 p 

(27T) 4 



Ag j (-P) 



/~2 _|_ ^2 _ — 777,2 ^2 _|_ ^2 _ 2&jj 



mi 



(76) 



The ie is assigned to fco by the Landau-prescription. Exploiting that p(p) ~ 5(p 2 — m 2 ) one can 
write 



S a b(k) 



dp 



+ 



where M 2 



fc 2 - 2fcp — M 2 A; 2 + 2kp + M 2 
6 — m 2 . Finally performing a ±fc/2 shift the result is 

d 4 p Ag? (p - fc/2) - Ai j (p + fc/2) 



(77) 



Sab(k) 



{2k)* 



2kp + M 2 



(78) 



Finally, introducing for the Fourier-transform of the left hand side of Eq. fl7l| ) the represen- 
tation 

W> = / ^AW(k,p), (79) 



we obtain from Eqs. fl7l| ) and (|7i 

2fcpA( 1 £ )(fc,p) = -Xa^ik) [AW(p + fc/2) - A® (p - fc/2) 
(2fcp + M 2 ) A^(fc,p) = -^#$i(fc) [Ai 0) (p + fc/2) - a2 } (p - fc/2)' 
which exactly coincides with Eq. (| 



(80) 



B Appendix 

An alternative approach for calculating characteristic quantities of real-time correlation func- 
tions is based on real-time dimensional reduction ]23| , 24] and solving the resulting classical 
effective theory relevant for the modes with high occupation numbers. On-shell as well as the 
off-shell damping rates were successfully reproduced [25, in the symmetric phase of 4> theory 
using this approach. In this appendix we will discuss the application of the classical approach 
to the 0{N) model and compare its results with the exact one-loop dynamics. The Lagrangian 
of the classical O(N) theory has the following form: 



The classical equation of motion corresponding to this Lagrangian is: 

(a 2 + m 2 (A))l> a + ^l» a (l> 2 )+j a = 0. 



(81) 



(82) 



The external currents j a (x) were introduced in order to prepare the derivation of the classical 
response theory. They should not be confused with the induced currents J a appearing in the 
effectice quantum equations of motion. The classical mass m c /(A) is different from the mass 
parameter m of the quantum theory (see Eq. (|3|)). The same is true for the coupling A c ;. The 
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classical mass parameter depends on the ultraviolet cutoff A. The results of the classical and 
quantum calculation could be matched by a suitable choice of m c i(A) and A c ;. In particular 
we shall see below (eq.(|9l|)), that the ultraviolet divergencies of the classical theory can be 

|| . Therefore we will 



eliminated if the divergent part of m c \ (A) is suitably chosen 1 23 , |24], |2!_ 
separate out the divergent part from the classical mass and write (A) = rrij. + 5m? (A) []. 
In the broken symmetry phase one separates out the condensate $, 



and the equations of motion read 



(d 2 +m 2 T + ^ 2 )$i + + ^1 + + 6m 2 $x + ji 

2 b 2 6 



{d 2 



A, 



6 ' * 6 v 0/ 3 
In addition the following initial conditions are imposed: 

$ a (t = 0,x) =F a (x), ^$ a (t,x) 



t=0 



P«(x). 



(83) 

(84) 
(85) 

(86) 



The expectation value of some quantity O (e.g. some correlation function ) is obtained by 
averaging over the initial conditions with the Boltzmann factor determined by the classical 
Hamiltonian H c i(P a , F a ,<&) corresponding to (|81~D 



< O >-- 



1 

z 



DF a DP a Oexp(-f3Hci(P a ,F a ,$)) 



(87) 



Z = / DF a DP a eM-PH cl (P a , F a , $)), 



The explicit form of H c i(P a , F a ,Q) is 



H d {P a ,F a ,$) 



d 3 x 



\PI + \^F a f + \{ml + ^)Fl + \{m 2 T + ^ 2 )F? + 



where we have separated out the classical condansate $ since no averaging over classical con- 
densate is understood. Since the classical condensate is separated both from the dynamical 
fields and from the initial conditions the following equation holds: 



< * a (t = 0,x) >=< Fi(x) >= 0. 



(90) 

Using the explicit form of H c i (Eq. (|89[)) this results at one-loop level in the following equation 
for the classical condensate <&: 



A 



(77^(A) + ^3> 2 )$ + ^$ 



A, 



6 



6 



d 3 p 

W 3 



1 



p 2 + m\ + ^-$ 2 



+ (N-1) 



1 



p 2 + nij, + 



0. (91) 



Choosing 5m? (A) to cancel the linearly divergent piece on the left hand side this relation reduces 
in the high temperature limit (mj>, A c /$ <C T) to 



m| + ^$ 2 
o 







(92) 



3 The divergent part of the classical mass parameter will be treated as interaction, similarly to quantum field theory 
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and 

Sm 2 = -^(N + 2)^. (93) 



Our procedure of solving the classical theory perturbatively closely follows Ref. [25]. Equations 



(84), (|85j ) are rewritten in form of the following integral equations: 

Qxfaj) = + J d 4 x'D R (x - x'X^i^a) + + ?f**i + Sm^+h) 

$i(x,j) = $?(x) + J d 4 x'D R (x - a0(^f$i(*2) + Y**i** + <5m2 **+i*), ( 94 ) 

where 

^-^) = -*w/^^ < 95 > 

is the classical retarded Green function p^] with w g = \J q 1 + m*, where raf = + 4^-<|> 2 = 
A^-^ 2 is the Higgs field mass and mf = + ^-$ 2 = (we have used Eq. (j92"D). Furthermore, 
<J?^ are the solutions of the free equations of motion and j in the argument of $ a refers to 
the functional dependence on j = Following Ref. one introduces linear response 

functions 



Using the integral equation (|94f) one can derive a coupled set of integral equations also for the 
linear response functions H R by functional differentiation of eq. (|94|) (see Ref. p5| for details). 
These integral equations can be solved iteratively in the weak coupling limit. When solving the 
equations, it is important to exploit the fact that only diagonal components of H R have terms 
O(X ). The classical retarded response function is the ensemble average of H R with respect the 
initial conditions : 

G&{x - x') =< H% b (x -x')>. (97) 
By eq.(p4|) one easily finds that it satisfies a Dyson-Schwinger equation of the general form 

') = D a R {x - x')5 ab + J d%d*y'D a R (x - y)5 ad U^ c (y - y')G? b {y' - x'), (98) 



G^h(x — x 



where n^ 6 (y — y') is the classical self-energy. 

In the perturbative expansion the average is done with the free Hamiltonian and therefore 
all thermal n-point functions are expressed as products of the two-point function of the free 
fields (solutions of the free equation of motion). The two point function of the free fields reads 



||, |3| as 

A^°(x - x') =< $°(x) a $°(x% >°= T f -^ e ^( x - x ')i_ cosK(t _ t ')) (99) 

J (27T) U q 

This classical two point function is analogous to the free two point function of the quantum 
theory A^J(x, x') Using eq. ( p4|) one gets the following self-energies for the Higgs and Goldstone 
fields at leading order in the coupling constant A c /. 

nf 1 Kk)=5m 2 + ^^A^(0,0) + 

{\ d <bf j dtd^ t -^[^ fi (y)D R {y) + E^l A f°(y)DUy)), 

18 



H^,k) = 5m 2 + -^5>f°(0,0) + 

a 

. - 2 

(^p | did 3 ye^- ik y(Af' Q ( 2/ ) J Dl,( 2 /) + Af' Q ( y )^( 2/ )), (100) 



Using the explicit form of A^'°(0, 0) and Eq. ( |93|) one can easily verify that all divergencies 
present in the above expression cancel. Furthermore, in the high temperature limit and at lead- 
ing order in the coupling constant only the last terms contribute in the expression of Ilf^uj, k) 
and H^(oj, k). One has to evaluate integrals of the following intrinsic form 

f d 3 p f dp dp' p a (p ,Lj a )p b (p' ,LU b ) d . cl , 

I dt l d Ve yD ^ A " {y)= J J^fJ u J -p -p' + i e {n {P ° ) + n (Po)) ' 

(101) 

where oj a = \/p 2 + m 2 and uj b = J (p + k) 2 + m 2 . The above expression coincides with the 



result of the quantum calculation (74), except the fact the there is no T = contribution and 
the Bose-Einstein factors are replaced by the classical distribution: T/pq. Then it is easy to 
write down the explicit expression for the classical self-energies using the formal analogy with 
the result of the one-loop quantum calculations. For example for the classical on-shell damping 
rate of the Goldstone modes one easily gets the following expression: 

cl _ Imn^ = |k|,k) = A 2 * 2 ff^\ d = Xl^T ( 4jkf \ 
1 W 2|k| 288vr 2 |k| 2 7j§ P W 288vr 2 |k| 2 V m l ) ' 

(102) 

Now let us discuss the correspondence between classical and quantum calculations. It was 
shown in Refs. [23, 24, 25, |2(| that the result of classical calculations can reproduce the high 



temperature limit of the corresponding quantum results if the parameters m 2 ; (A) and A c / are 
fixed to the values determined by static dimensional reduction. For our case this implies: 

A /T 2 ATX 

rn 2 cl (A)=m 2 + ^(N + 2) —-—,), X d = A, (103) 



6 v ' V 12 2?r 2 J ' 

where m and A are the mass and the coupling constant of the corresponding quantum theory. 
From the above equations it is easy to see that the divergent part in m d coincides with 5m 2 
determined from the ultraviolet finitness of the classical result (cf. Eqs. (|9ll) , (|93|) and ( |100f ) ) 
and mi is the high temperature limit of M\ (see Eqs. (92) and (|l3|)). For small values of |k| 



(|k| <C mi) the logarithm in this expression can be expanded and one obtains the result of 
Eq. (|40|), which fails to reproduce the result of the quantum calculation. The same is true 
for the whole imaginary part of the Goldstone self-energy. The high temperature limit of the 
imaginary part of the Higgs self-energy (see Eq. ([?(])), on the other hand, is well reproduced 
by the classical theory. 

This result can be easily understood by looking at the effective quantum equation of motion 
(|i~T|). in the equation for the Goldstone fields the induced current Jj has a non-local contribution 



from loop momenta p ~ M 2 /|k| >> T (cf. (32)). However, no such term is present in the 



corresponding classical equation of motion on one hand and the classical theory cannot describe 
fluctuation with wave length much smaller than T~ l on the other hand. The induced current 
J\ in the effective equation of motion for the Higgs fields receives non-local contribution only 
from the loop momenta around p ~ M\ which can be described in the framework of the classical 
theory. 
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C Appendix 



In this Appendix we illustrate the steps of evaluation of the integrals in Eq. (29). 
Imaginary parts. As example of the evaluation of a relevant integral we discuss in detail 

i-niw*, m) = im / ^ A(»»( y - t /2)-A<°)(p +t/ 2) i (104) 

where the propagators can be either of type 11 or ii In order to implement the Landau- 
prescription we transform away any k dependence from the propagators by shifting the p integral 
by ±/c/2. Then we use 

lim Im = — ne(a)8(x), (105) 

a^o x + ia 

and finally shift the integrals back. We write the 4D integration measure as 



d 4 p 1 



oo oo 1 



dpp 2 I dpo I dx, (106) 



(2vr) 4 8vr 3 



where x = pk stands for the cosine of the angle between the spatial momenta. The x integration 
is trivial because it appears in the Dirac-delta arising from the application of the principal value 
theorem 

l 

J dx5(p k -p\k\x) = 4r©(p|k| - \p k \). (107) 

Using the explicit form of the propagators (see Eq. (|25|)) and the identity 0(a>) + n(|cj|) = 
e(w)(l + n(u>)) we find 



CO 



-hnR x {k,M) = —j - JdppJdp e(po-^-)e(p + f)5(p 2 -S 



™(Po-y) -"(Po+y, 



(108) 

where a = p|k|/|fco| an d S 2 = p 2 + M 2 — k 2 /A. The po integration over the Dirac-delta gives a 
constraint for the p integration of the form 



S < a => p -J < M . (109) 

k 4 

This can be fulfilled only for k 2 > AM 2 (above the two-particle threshold), or for k 2 < (Landau 
damping). After elementary algebra one can establish the value of the sign functions and one 
arrives at the formula appearing in Eq. (p3|). 

A similar analysis can be performed for Imi?j, however, in this case it proved to be more 
convenient to start from the equivalent form 

After implementing carefully the Landau prescription and performing the x integration as de- 
scribed before, we arrive at 

oo £>+ 

-lmRi(k,M) = — — J dppj dp e(p ) e(p -k ) 5(pl-p 2 ) [n(p -k ) - n(p )] , (111) 
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where 



b± 



'-M 2 |k| 
± p — . 

2k k 



112) 



The po integration over the mass-shell delta-function again restricts the domain of integration 
in the p integral: 6_ < ±p < 5 + , which, however, does not restrict the possible values for k 2 . 
After the solution of these linear inequalities and the analysis of the sign functions we arrive at 
the result appearing in Eq. (|33|). 



Real parts. The relevant integrals in Rei?i are 

d 4 p A<°>(p) 



Re 



(2vr) 4 pk±k 2 /2' 



(113) 



where the propagator can be either of type 11 or ii. With their help we find Rei?i = I + — I . 
The real part comes from the principal value integration. We decompose the integration measure 



as in the calculation of the imaginary parts. When we use Eq. (25) for the propagator, the value 
of po is fixed by the delta function. The x integration can be performed as 



dxV 



2r>M) — 2p|k|x ± k 2 p|k 



1 ,w 2p|k| 
arth(- 



'2poko ± k 2 



(114) 



where arth(rr) = 1/2 In |(1 + x)/{l — x)\. Using the properties of the absolute value we find 

7 + (Ad,k) = -J-(-fc ,k). (115) 



Then we directly arrive at Eq. (H). 
We write Rei?,- in the form 



ReRi = Re 



Ag J (p) 



d 4 p 

(2vr) 4 2pk + k 2 + M 2 



Re 



d 4 p 

(2vr) 4 2pk -k 2 + M 2 ' 



(116) 



Then the previous scheme of calculation goes through directly. 

The only problem still to be discussed is the zero temperature contribution, which diverges 
logarithmically. The regularization and renormalization prescriptions are better formalized in 
the language of the propagators and couplings. Using the results of Appendix [A| to go over to 
the perturbation theory, we have to evaluate 



R 



d 4 p 

WT 4 



G c (p)G u (p - k) - G < {p)G < {p - k) 



iC 



(117) 



The first term at T = is the usual time ordered product; with finite four-dimensional cutoff 
one finds 



d 4 p 



G c (p)G v {p-k) 



1 



16vH 



i 

1 + J dx In 





\k 2 x{\ — x) — m 2 x — m^l — x) 
A 2 



;ii8) 



The divergence is canceled by the coupling constant counterterm. In modified minimal subtrac- 
tion (MS) scheme it reads 



16vr 2 



1 , , \k 2 x(l — x) — m 2 x — mo(l — x) 

dx In- i — 

fi 2 



(119) 
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In the second term of Eq. (117) at T = we can use G < (p) = @(—pq)(2tt)5(p' 2 — m 2 ). 
Because of the delta functions this piece yields finite result (as is expected by the arguments of 
the renormalizability) 



g^Ip-*). "L {£-¥ - e (f -^").» »i = ™ 2 = ™ (120) 

( 27 87r (l-^-)6(A: 2 -m 2 ), if m x = 0, m 2 



m. 
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